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J_h . Unitarity corrections to the BFKL description of high energy hard scattering are viewed in 

£L|- large N c QCD in light -cone quantization. In a center of mass frame unitarity corrections to high 

energy hard scattering are manifestly perturbatively calculable and unrelated to questions of parton 
saturation. In a frame where one of the hadrons is initially at rest unitarity corrections are related 
to parton saturation effects and involve potential strengths ~ 1/g. In such a frame we describe 
the high energy scattering in terms of the expectation value of a Wilson loop. The large potentials 
A^ ~ 1/g are shown to be pure gauge terms allowing perturbation theory to again describe unitarity 
^ ' corrections and parton saturation effects. Genuine nonperturbative effects only come in at energies 

, well beyond those energies where unitarity constraints first become important. 

I. INTRODUCTION 

o 

Q\ ■ Our object in this paper is to describe a one (hard) scale high energy scattering in QCD at those energies where 
unitarity corrections begin to be important. In particular, we focus on the question of the perturbative calculability 
of such processes. As our prototype reaction we consider high energy heavy onium-heavy onium scattering where the 
hard scale is given by the inverse onium radius, the only relevant transverse momentum scale for the total cross section 
Oh, or for a near forward elastic scattering. We could equally well focus on more directly physical processes such as the 
total virtual photon- virtual photon cross section jj]|2| 7*(Q) + 7*(Q) — * hadrons, or on associated jet measurements 
either in electron-proton collisions J^^] or in proton-antiproton collisions ||, however unitarity issues are especially 
^ ■ sharp in onium-onium scattering since the onium is a bound state in QCD. 

km | The energy dependence of an onium-onium collision is governed by the BFKL equation , appropriate for single 
scale high energy collisions. If the onium is sufficiently small one can neglect diffusion and running coupling effects 
and use the BFKL equation up to those energies where unitarity corrections become important. We shall also use the 
large iV c limit of QCD where one can view the light-cone wavefunction of a high energy onium as a collection of color 
dipoles |i"o|-|l2]]. An onium-onium collision, in the BFKL approximation, consists of a scattering of a single dipole in 
each of the onia by means of one gluon exchange. 

The BFKL approximation gives an S'-matrix for onium-onium elastic scattering at impact parameter b and rapidity 
Y behaving as 
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a behavior which violates the unitarity bound \S(Y, b)\ < I when Y becomes very l arg e. The picture of unitarity 
corrections is very different in different Lorentz frames in light-cone quantized QCD |jl^] in light-cone gauge. In the 
center of mass frame unitarity corrections become large when dipole densities, and corresponding field strengths A^, 
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are relatively small. Unitarity bounds are guaranteed by multiple dipole-dipole scattering [|10|,|l4| in much the same 
way as occurs for nucleus- nucleus scattering in the Glauber model. One difference, as emphasized below, is that 
the multiple dipole scattering must be done for a fixed wavefunction configuration with a sum over the different 
configurations performed as a last step. The higher number of scatterings correspond to higher orders in l/N c . In the 
center of mass system unitarity corrections are quite straightforward and numerical simulations have already been 
done jL4|. Unitarity limits are enforced at rapidities well below those where field strengths become large and parton 
saturation effects are important. 

It is instructive to view unitarity corrections in a different Lorentz frame |jl3f where one of the onia, say the 
right-mover, carries almost all the available rapidity with the left-moving onium having a high enough rapidity to be 
relativistic but not so much that higher (gluonic) left-moving Fock components must be considered in the scattering. 
In this frame, called the L-frame, the scattering is of a right-moving onium with a fully developed wavefunction of 
many gluons colliding with a left-moving onium whose wavefunction is simply a heavy quark- antiquark pair. The 
S-matrix can be written as an average of Wilson loops in the right-moving state as indicated in (30). From this 
representation it becomes clear that unitarity corrections become important in the L-frame when ~ 1/g. In this 
frame unitarity is associated with high potential strengths and significant parton saturation. However, it is suggested 
that these large values of A^ do not yet signal nonperturbative QCD effects since the large terms in are purely 
gauge terms. In order to see genuine nonperturbative QCD effects one must go to rapidities well beyond (about twice 
as large) those rapidities where unitarity corrections first become important. 



II. UNITARITY CORRECTIONS IN THE CENTER OF MASS FRAME 

In this section, we review high energy heavy onium - heavy onium ]Toj — |T^| scattering in large N c QCD beginning 
at moderate energies where the two gluon exchange approximation is valid, through the energy range where the 
BFKL equation governs the scattering, up to energies where unitarity corrections to the BFKL equation become 
important, into the energy domain where unitarity corrections are very strong and the forward S-matrix elements 
become nearly zero (blackness), and ending at energies so high that the whole perturbative picture breaks down and 
high field strengths (parton saturation) occur. 

Suppose <&(a;oi, z) is the square of the heavy quark light-cone wavefunction of the onium with xo± = — x | the 
transverse coordinate separation between the heavy quark, at x , and the heavy antiquark, at Xj and with z the 
fraction of the onium's longitudinal momentum carried by the quark. Then in the two gluon exchange approximation 
the forward elastic high energy scattering amplitude is 



A=- J d 2 x 01 d 2 x' 01 J dzdz'$(xoi,z)$(x' ol7 z')<T DD (xoi 7 x' Q1 ) (1) 
where the dipole-dipole cross section, ctdd, is 

o-dd(xqi,x' 01 ) = 2ira 2 x 2 < [l + ln(x > /x < )] (2) 
with x < the smaller of xqi and x' 01 and with a;> the larger of these quantities. The onium-onium cross section is 

o- = 2ImA (3) 
which is independent of energy at high energies. The a in (2) is a(R) where 



R = - j d 2 x 01 J dzx 01 $(x i,z) (4) 

is the average light-cone radius of the heavy quark part of the onium wavefunction. If RAqcd is sufficiently small, 
running coupling effects will not be important in onium-onium scattering up to energies where saturation sets in and 
the whole perturbative picture of high energy scattering breaks down. 

If Y = ln(s/M 2 ) is the relative rapidity of the two colliding onia then the two gluon approximation breaks down 
when aY is of order 1. In this regime the BFKL equation, an equation which sums all (aY) n terms, governs high 
energy onium-onium scattering. In the dipole picture of high energy scattering one writes the scattering amplitude 
at impact parameter b as jT^] 
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A(b, Y) = — i I d xo\d x' Q1 j dzdz'$(xoi,z)$(x' 01 ,z')F(xQi,x' 01 ,b,Y) (5) 
where in the BFKL approximation F = and 



-i />00 7 J I 

F (1) = -^ — "4 d 2 bm(x 01 ,x, h, Y/2)n(x' 01 , x',\b- h^Y /2)<j dd {x, x') 

Z Jq X X 



(6) 



n(xoi, x, bi, Y/2) is the number density of dipoles of size x whose center is a distance bi from the center of the heavy 
quark- heavy antiquark system, xqi. The forward elastic amplitude is 



A(Y)= J d 2 bA{b,Y). (7) 

For x/Moi/&< 1 011,0 

, xqi /16b 2 \ exp{(gp - 1)Y - a(Y)/2 In 2 (16b 2 / x 01 x)} 

with a P - 1 = (4:aN c /n) In 2 and a _1 (^) = 7ajv ^ (3)F , Using (2) and (8) in (6) gives 

r nn r' / ^(\h 2 \ expUap - 1)Y - fln 2 (^^-)} 

when xoi/b,x' 01 /b <C 1. The resulting onium-onium cross section 

e (a P -l)Y 

a = l&nFfa 2 —^ (10) 
^laN c ((3)Y 



comes from impact parameters satisfying In y-^z) ~ y/14aN c £(3)Y/iv. From (3), (5), (7) and (9) one obtains 

rf 2 6 — ' V— [I a iV c C(3)y] 3 / 2 & 2 

so long as R 2 /b 2 <C 1. Since 



J^ = l-ReS(6), (12) 

with S(b) the 5*-matrix for onium-onium scattering at impact parameter b, (10) and (11) violate the unitarity bound 
4%r < 2 when aY" becomes large. 

The dipole picture of BFKL behavior is very simple. At high energies each of the onia consists of a heavy quark- 
antiquark pair along with a (generally) large number of soft gluons ordered in longitudinal momentum. In the large 
N c limit each gluon can be viewed as a quark- antiquark pair and the onium wavefunction then can be viewed as a 
collection of color dipoles. In the scattering of two onia, in the BFKL approximation, a single right-moving dipole 
from the right-moving onium scatters off a single left-moving dipole in the left-moving onium, with the two gluon 
exchange cross section given in (2), leading to the parton-like expression (6). 

Roughly speaking, one expects unitarity corrections to the BFKL approximation to become important when -fc 1. 
From (2) and (6) this happens at values of Y when n ~ 1/a and one might expect that the leading logarithmic 



3 



approximation to the onium wavefunction would start to break down. However, it is straightforward to see that when 

da -i 

x 2 n~^-. (13) 
ab 

Since x/b -C 1 for those dipole sizes dominating the integral in (6) the onia wavefunctions are still quite dilute when 
unitarity corrections begin to become important. The situation is somewhat analogous to the scattering of two large, 
but dilute, "nuclei." If the nuclei have atomic number A and radius R then, in the single scattering approximation 

ctaa = A 2 a = 4kR 2 P 2 ^J) , (14) 
where <jq is the nucleon-nucleon cross section and p is the packing fraction defined by 



Act, 



3/2 



n 



(4/3)7Ti? 3 ' 



(15) 



From (14) one sees that the unitarity limit, AttR 2 , is reached for cjaa for very small packing fraction if R 2 /(Tq 1. 
In this nuclear example it is clear how to impose unitarity. Since the nucleus is dilute one need only do multiple 
scattering quantum mechanically to obtain a correct forward scattering amplitude. 

The unitarity limit is reached for onium-onium scattering in the center of mass frame also when the dipole densities 
in the onia are small so that the two dipoles involved in the scattering can be viewed as not interacting with any other 
dipoles. It would seem that unitarity can be imposed simply by including two and more scatterings of left-moving 
dipoles with a similar number of right-moving dipoles. And this is indeed the case 0,|l| . There are, however, two 
differences from the "nuclear" example, (i) The dipoles have various sizes whereas the nucleons in our nuclear example 
all have the same size. This is a technical problem which makes an analytical discussion of the problem difficult but 
does not present any problems of principle pp| |. (ii) The wavefunction of an onium has large fluctuations in the 
number of dipoles present. Thus the single scattering term in the S-matrix is dominated by very different parts of 
the wavefunction than is the double scattering term which in turn samples a very different part of the wavefunction 
than the triple scattering term, etc. This causes the normal multiple scattering series (multiple pomeron exchange) 
to diverge factorially JTyJl^] . This means that one must sum the multiple scatterings of the dipoles in the colliding 
onia for a fixed wavefunction configuration and only after this sum has been done can the average over the different 
configurations of the onia wavefunctions be done. Schematically, 

s(Y,b) = Y / ^ fW p^y) as) 

where <j> labels configuration of dipoles in each of the colliding onia. That is <f> labels the number of dipoles in each 
onium along with the size of each of the dipoles and the impact parameter of the center of each of the dipoles with 
respect to the center of the heavy quark-antiquark pair. gives the probability of configuration <f> while /(</>) is 
proportional to the probability of a dipole-dipole scattering in the configuration <f>. b is the impact parameter of the 
onium-onium collision. If one expands the exponential in (16) and interchanges the sum with tf> 

oo 

S(Y,b) = J2(-l) n S n (Y,b) (17) 

n=l 

with 

gives the usual multiple scattering series. However, simple models |10 and numerical calculations in QCD Q show 
that S n ~ n\ so that the multiple scattering series does not converge. The problem is that rare configurations which 
have a large value of f{<j>) contribute very much to S n , for large n, but very little to S(Y,b) as given in (16). 
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Numerical simulations have been carried out within the dipole framework of high energy onium-onium scattering 
which follow the BFKL behavior given in (10) and (11) into the regime of ^-values where unitarity corrections due 
to multiple scattering are important. These simulations suggest that Q 

S(Y,b) ~ e'< Y ' Y ^ 2 (19) 

for Y much greater than Y(b) 7 the value of Y at which unitarity corrections become important. The behavior given 
in (19) is a much slower decrease of S with Y than would be expected from an eikonal picture and corresponds to the 
S matrix being dominated by rare configurations consisting of many fewer dipoles than the average at the Y value 
being considered. The behavior (19) is also expected on theoretical grounds. 

If Yq is the rapidity at which unitarity corrections begin to be important for b = collisions then at Y « 2Yq one 
expects perturbation theory to break down completely and a new regime of high field strength QCD to emerge. This 
can be seen from the following rough argument: In (6) neglect the integrations over dxdx 1 and set d 2 b\ = R 2 for a 
b = collision. (The neglected integrations only give logarithmic prefactors in any case.) From (6) it is clear that 
unitarity corrections become strong when R 2 n(Yo/2) ~ 1/a. However, at this value of Yq the onium wavefunction is 
still relatively dilute. Indeed, if one scatters a single low momentum left-moving dipole on a fully developed right- 
moving onium state of rapidity Yq/2 the cross section is of order a. The interaction of the fully developed right-moving 
onium state with a slightly left-moving dipole is the same as the interaction with a slightly right-moving dipole which 
could be part of the onium wavefunction. It is only when R 2 n f» 1/a 2 that a given low momentum dipole in an 
onium wavefunction begins to strongly interact with some of the dipoles of higher rapidity. Thus in a center of mass 
collision at Y = 2Yq all the soft dipoles in both the right-moving and left-moving onia suffer a strong interaction and 
a good portion of these dipoles will appear as freed gluons in the collision leading to a high field strength initial state, 
R 2 F fll , <~ 1/g, of the collision. However, at such high field strengths our control over the collision is lost because 
perturbation theory breaks down. (This regime is usually referred to as the saturation regime because a limit to the 
growth of parton densities is expected to occur.) Whether there is some classical, or semiclassical, solution which 
governs such collisions is an interesting and challenging problem which, however, goes far beyond the purposes of the 
present work. 

III. SCATTERING IN AN ALMOST REST SYSTEM 

In this section we shall view the forward scattering of two onia in a frame where one of the onia, say the right- 
moving onium, has almost all the available rapidity while the left- moving onium has only enough rapidity to make 
it move relativistically, but not enough so that one is required to add gluons to its wavefunction p3| . Thus the 
left-moving onium in this almost rest system, which we shall call the L-frame, consists only of a heavy quark and a 
heavy antiquark, a single dipole. 

A. The BFKL Approximation in the L-Frame 

The BFKL approximation is straightforward to implement in the L-frame. In (6) one just makes the replacement 

n(x' 01> x', |b - fail, Y/2) - n(x' 01> x', |b - b x |, Y) (20) 



n(x i,x,b 1 ,Y/2) -> xS(x - SGoi)5(bi)- (21) 

Eq.(21) simply says that the left-moving onium consists of a single dipole, the heavy quark-heavy antiquark pair. The 
result (9) emerges as is easy to check. 

However, it is useful to view the process in a somewhat different way [^6|-p^|. Let us view the interaction of the two 
onia as they pass each other as due to the interaction of the left- moving heavy quark-antiquark pair with the color 
field of the right-moving collection of color dipoles making up the right-moving onium. Thus as the onia pass each 
other we consider a fixed configuration of dipoles in the right-moving onia. Each dipole gives a classical color field 
(We take traces over the colors in the individual dipoles later.), and because of the large N c approximation the color 
fields coming from the different dipoles do not interfere. The interaction with the left-moving heavy quark-antiquark 
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pair with a dipole in the right-moving onium of size x' at impact parameter bi from the center of the left-moving 
heavy quark- antiquark pair is 



1 



trP exp ig <z> A ■ dx 



(22) 



where the integration goes around a closed Wilson loop in the x±,x^ plane as illustrated in Fig.l. Indeed, explicit 



evaluation identifies W 2 with — j^f(b 1 ,x',x 01 ) of Ref. 



1$. The in (22) is the field due to the color dipole of size 



x' and the subscript 2 in (22) means that we are only taking the order g 2 term. The exact values of x_ and x r _ for 
the Wilson loop do not matter so long as x r _ — x_ > R and x r _ > 0, x e _ < with neither x r _ or x e _ being too close 
to zero. (We suppose the centers of the two onia pass each other at x+ = x_ = 0.) The field, not yet counting color 
factors, due to the dipole x' , consisting of a "quark" at x[ a and an "antiquark" at x[ bl is [Rl| 



S(x-) ( {x 



M x )=9^^(^^)^ (23) 



\ 2 



^ 4tt ^\{x~x' a ) 2 

in covariant gauge and 



M x ) 



47T 



e(x- 



( x x 'a)ji ( x X 'b)ji 
(x-xj,) 2 " (X-X'J2 



(24) 



in principal-value light-cone gauge. In (23) rj is such that rj ■ v = i>_ for any four-vector i> M while xj^ in (24) means 
that only the fi = 1, 2 components are nonzero. From the forms of A^ in (23) and (24) it is clear that the exact values 
of x r _ and x e _ will not matter when evaluating the Wilson loop in (22). 
When xoi/b, x' 01 /b 1 it is straightforward to check that 



F^(x 01 ,x' 01 ,b,Y) = J ^^d 2 b 1 n(x' 01 ,x',\b-b 1 \,Y)W 2 ( } c',x 01 ,b 1 ) (25) 

agrees with (9). Eq.(25) is in fact true for all b in the leading logarithmic (BFKL) approximation. (Also, the lowest 
order approximation, (21), for n allows one to identify W 2 as the dipole-dipole scattering amplitude at a definite 
impact parameter.) One can use (25) to find out what are the typical values of the fields A^ in the right-moving 
onium. From (22) and (25) 



^trUgjA-dx^j l=- 2 J d2x '°iJ Q dz'$(x' 01 ,z')FW{x' 01 ,x 01 ,b,Y) (26) 

where the left-hand side of (26) is the expectation in the right-moving onium state of the square of the loop integral. 
This is now the field due to all of the dipoles in the right-moving onium. In the BFKL approximation (26) is a gauge 
invariant equation. Using (9) one finds 




A . \ = W <!XP{( V - . ,27) 



In light-cone gauge the integration of the A^ field goes over distances xoi so that the integration picks out those 
parts of A^ having wavelength xqi or greater. Define the left-hand side of (27) to be -^-x\ x g 2 ((A^) 2 ) . (The factor 
of 4 is included because, according to (24) there are equal contributions coming from x_ < and from a;_ > in the 
integral on the left-hand side of (27).) Then, 



ln 



/ 16b 2 \ CX P 
^Rx i) 



\a P -l)y- § ln 2 (^l) 
[|aiY c C(3)F] 3 / 2 



(28) 



6 



gives the average value for the square of a fixed color component and a fixed transverse spatial component of 
in light-cone gauge when the field is measured a distance b from the center of the fast moving onium and if the 
measurement averages over transverse distances of size Xoi- When aY is of order 1 the field is of size g. Taking 
xqi — 2R and using (11) 



so that the unitarity limit, in the BFKL approximation, is reached when typical field strengths reach A <~ {gy/Nc) -1 ■ 
Hence, when light-cone field strengths reach values of the order of l/g one expects strong unitarity corrections to 
occur. Thus, in the L-frame the onset of unitarity corrections to the BFKL approximation occurs at exactly the same 
rapidity that field strengths reach size l/g. This is in contrast to the center of mass frame where unitarity corrections 
become important when field strengths are of size g°. 



B. Unitarity Corrections in the L- Frame 



When field strengths reach a size l/g the complications of doing calculations in the L- Frame are of two varieties. 
First, one can no longer expect the two gluon approximation to be valid for connecting the right-moving system of 
a heavy quark-antiquark pair and its accompanying right-moving gluons to the left-moving heavy quark-antiquark 
pair. Secondly, one can no longer expect the right-moving system to consist only of dipoles. When A 1 ^ ~ (gy/Nc) -1 
the large N c expansion begins to break down and one can expect more complicated color structures, quadrupoles, 
sextupoles, etc. to become important in addition to the color dipoles corresponding to the BFKL approximation. In 
the center of mass system the breakdown of the large N c expansion shows up in the multiple scattering in (16) and, 
more explicitly, in (17) and (18). 

We continue to find it convenient to view the scattering as the interaction of the left- moving heavy quark-antiquark 
pair (dipole) with the color field of the right-moving system. Thus from (5), with S(Y, b) = 1 + iA(Y, b) 



S{Y,b) = ±r j d 2 x al j\ Zl <Z>(x al , Zl )tr(Pe l9 $ A ^) b (30) 

where the contour of integration for dx^ is the same as before, and shown in Fig.l, while the average in (30) means 
that A^ is the field in the high momentum right-moving onium at fixed impact parameter b from the center of 
the right-moving heavy quark- antiquark pair and with averages being done over all the numbers and distributions 
of gluons (dipoles, quadrupoles, etc.) of the right-moving system making up the wavefunction of the onium. That 
the Wilson loop correctly expresses the scattering amplitude is apparent in covariant gauge where the interaction 
with the left-moving quark-antiquark pair is causal and occurs with the horizontal lines shown in Fig. 1. Since the 
left-moving heavy quark-antiquark pair is frozen, in transverse coordinates, during the passage through the field due 
to the right-moving system, the integral in (30) just represents the eikonal interaction of that field with the heavy 
quark at x and with the heavy antiquark at x^. 

In the BFKL approximation the fields from the right-moving system come from independent dipoles and (30) is 
replaced by (25). In general (30) shows that when one reaches rapidities where A 1 ^ ~ (g^/Nc)^ 1 , and unitarity 
corrections become strong, what is being determined in the scattering is the expectation of a Wilson loop. The 
expectations of the various sized Wilson loops contain the essential information of QCD. Very high energy scattering 
allows one to test QCD in a regime where g is small, corresponding to the small onium size, but where gA is not 
small. 

Suppose we take Y = Yo, the rapidity at which unitarity corrections begin to be important for a zero impact 
parameter onium-onium collision. At such rapidities gA^ is of order 1 in (30), and one might expect that essential 
nonperturbative physics would be necessary to evaluate A. However, at a corresponding rapidity we found that 
perturbation theory works well in the center of mass system where gA^ is of order g. How does one reconcile this? 
The resolution is that the field A^ ~ l/g in (30) when unitarity corrections become important, but that the large 
part of this field is purely gauge. The potential A^ comes essentially from a superposition of independent fields of the 
type shown in (24). The potential given in (24) leads to a field strength 



2tt 



4) 2 . 



(31) 
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In reality the 5(x-) should be smeared out (T^JT^] over an extent Ax_ oc Re~ Vab where y a b is the smallest of the 
rapidities of the gluons making up the "quark a" and the "antiquark b" which is the source of (24) and (31). Thus, 
although the potential is of size 1/g when Y = Yq the F + i's which come from the different dipoles which serve 
as sources for A have differing longitudinal extent Ax_ so that there is no fixed coherence distance Ax- over which 
F + i is of size 1/g. This means that the 1/g part of is purely a gauge. (The largest number of sources having 
a fixed coherence length are those dipoles which are (partly) made up of gluons having rapidity on the order of 1. 
When Y = Yq there is on the order of iVo ~ 1/a such sources. These sources lead to a non-gauge part of A^ of size 
Ap, ~ s/Nq ■ g ~ 0(1) which is the maximum non-gauge field present at V = Yq.) Nevertheless, the determination of 
the Afj, to use in (30) does involve nonlinear (gauge) interactions because gA is not small. It is not a trivial problem 
to evaluate A M for an onium state having rapidity Y>Y . For the somewhat simpler problem of a large nucleus where 
A ~ 1/g comes from different (non-overlapping) nucleons, the corresponding pure gauge field, which could be used in 
(30) to give hadron nucleus scattering, has been calculated in Refs. [jljj and Jig] . 

Thus at those rapidities where unitarity corrections begin to be important the potential A M becomes of size 1/g 
but this large part of the potential is purely gauge. As one further increases Y the rate of growth of A M will be 
considerably slower than for Y < Yq in keeping with (30) leading to (19). This slowing down of A^ corresponds 
to "parton saturation" and in principle is perturbatively calculable. In the L-frame parton saturation is intimately 
related to the rate at which unitarity (blackness) sets in as Y becomes greater than Yq juj). However, as we have seen 
in Sec. 2 one can calculate unitarity corrections in the center of mass system where field strengths remain weak and 
saturation effects are not important. 

Finally, we have argued that at values of Y on the order of Yq, where unitarity corrections first become important, 
the wavefunction of a heavy onium is still a perturbative object even though A^ ~ 1/g. Thus the general probe 
described by the Wilson loop in (30) and corresponding to scattering by an undeveloped onium state does not require 
knowledge of nonperturbative QCD. 

However, suppose we scatter, at zero impact parameter, a right-moving onium, having Y « Yq, on a left-moving 
onium, having Y m Yq. In the central unit of rapidity there are on the order of 1/a gluons, corresponding to 
A ~ 1/g, which will be freed during the collision. These gluons lead to field strengths F^ u ~ 1/g and now genuine 
nonperturbative QCD effects can be expected. Thus the genuine nonperturbative effects appear to be more connected 
with production in the collision process than with the wavefunction of a high energy hadron, or at least they show up 
at lower rapidities in the collision than in wavefunctions entering the collision. 
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